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ABSTRACT 

The Li abundances of the two components of the very metal-poor ([Fe/H]= —2.5) 
double-lined spectroscopic binary G 166-45 (BD+26°2606) are determined separately 
based on high resolution spectra obtained with the Subaru Telescope High Dispersion 
Spectrograph and its image sheer. Prom the photometric colors and the mass ratio 
the effective temperatures of the primary and secondary components are estimated 
to be 6350ibl00 K and 5830ibl70 K, respectively. The Li abundance of the primary 
(A(Li)= 2.23) agrees well with the Spite plateau value, while that of the secondary is 
slightly lower {A{L\)= 2.11). Such a discrepancy of the Li abundances between the two 
components is previously found in the extremely metal-poor, double-lined spectroscopic 
binary CS 22876-032, however, the discrepancy in G 166-45 is much smaller. The results 
agree with the trends found for Li abundance as a function of effective temperature (and 
of stellar mass) of main-sequence stars with —3.0 <[Fe/H]< —2.0, suggesting that the 
depletion of Li at reg~ 5800 K is not particularly large in this metallicity range. The 
significant Li depletion found in CS 22876-032B is a phenomenon only found in the 
lowest metallicity range ([Fe/H]< —3). 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances 
— stars: Population II — stars:individual(G 166-45) 
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Introduction 



Lithium is an important element as a diagnostic of the structure and evolution of low-mass 
stars. Since Li is easily destroyed by nuclear fusion in the hot interiors of stars where the tempera- 
ture is above 2.5 xlO^ K, stellar internal processes can be examined from observed Li abundances 
at stellar surface. Li depletion processes also play an important role in reconciling the serious dis- 
crepancy between the observed constant Li abu ndance, the so- c alled S pite plateau, of metal-poor 



turnoff stars (^(Li) = log(Li/H)-M2 ~ 2.2: e.g-. ISpite fc Spitd (|l982l ^: lAsplund et al.l (|2006l ^1 and 



the primordial Li abundance predicted from standard Big Bang nucleosynthesis (SBBN) models 
adopting the baryon density determined from observations of the cosmic microwave background 
(CMB) by the WMAP satellite ( A(U) = 2.72: e.g. Cyburt et al. 2008 ) . The discrepancy se ems 



to be larger at lower metallicity (jBonifacio et al.l 120071 : lAoki et al.l l2009l : ISbordone et al.l |2010| ) . If 
the stellar intrinsic Li depletion is in operation in the atmospheres of metal-poor turnoff stars, the 
observed Li abundance is lower than the initial value, which may have been closer to the CMB and 
SBBN prediction. 

Although no Li depletion is inferred for metal-poor stars on the main-sequence-turn-off by 
standard stellar evolution theories because of t heir shallow convectio n zones, evolutionary models 

redict significant decline 
additionally introduced 



that include effects of atomic diffusion (e.g., iRichard et al 



2002a 



Richard et al 



aub|) p 
Eo3) 



of surface Li abundances for metal-poor dwarfs , 

parametrized effects of turbulent mixing to their evolution model for [Fe/H] = —2.31 E|, t hough 
the understand i ng of the physics of the turbulent mixing is still insufficient. iKorn et al.l (|2006l ) 
and lKorn et al.l (j2007l ) investigated effects of these processes for the globular cluster NGC 6397 by 
examining small variations of Mg, Ca, Ti and Fe abundances in main-sequence turnoff, subgiant 
and giant stars. They suggested diffusion process during main-sequence phase and attributed the 
slightly lower Li abundances in turn-off stars than in subgiants to these effects. 

The most convincing evidence of stellar Li depletion in metal-poor dwarfs found so far is the 
large difference between Li abundances of the primary and secondary dwarf stars in the extremely 
met al-poor ([Fe/Hl= —3.6) double-li ned spectroscopic binary CS 22876-032, which is about 0.4 



dex (jGonzalez Hernandez et al.l 120081 ). Since the two stars should have shared the same chemical 
compositions at their birth, it is reasonable to presume that depletion processes have made the Li 
abundance difference during the stellar evolution. Such depletion is not expected for a star with 
the effective temperature of the secondary (Tefj = 5900 K) from previous studies of less metal-poor 
stars. Given that the mass of the secondary star is about 0.9 times that of th e primary, Li depletio n 
may operate m ass-dependently , which is consistent with the predictions of iRichard et al.l (|2005l ). 
as discussed by Melendez et al. ( 201ol ). 



The best approach to examine the hypothesis of the mass-dependent lithium d epletion is accu 



rate determinations of Li abundances for several metal-poor binary stars. Although iMelendez et al 



i[A/B] = log(iVA/iVB) - log(iVA/iVB)0, and logeA = log(iVA/iVH) + 12 for elements A and B. 
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([20101) suggested a correlation between Li abundance and stellar mass for metal-poor single stars, 
investigation of single stars which have different ages and distances would suffer from uncertainties 
in mass estimates. 

In this Letter, we report our chemical abundance analysis, in particular the Li abundances, 
of the two components of the very metal-poor, double-lin ed spectroscopic binary G 166-45. This 
object is selected from the sample of lGoldberg et al.l (j2002l ) who investiga ted the orbital para. meters 
of 34 binary stars found by the Carney-Latham proper motion sample (ICarnev et al.lll994l ). The 
chemical abundances of individual components of this object have not yet been investigated. We 
discuss in particular the implication of the small difference of the Li abundance between the two 
components found from the present analysis of high resolution spectra. 



Observations and measurements 



A high-resolution spectr um of G 166-45 wa s obtained with the Subaru Telescope High Dis- 
persion Spectrograph (HDS; iNoguchi et al.ll2002l ) on July 22, 2011. We applied the image sheer 
recently installe d in the spectrogra ph, which provides the resolving power of ii = 110,000 with 
high efficiency ( Tajitsu et al.ll2012l ). The wavelength range from 5100 to 7800 A is covered by 
the standard setup of StdRa. The exposure time is 60 minutes, resulting in the S/N ratio (per 
0.9 km s-i pixel) of 350 at 6700 A. 



According to the orbital parameters of the binary determined by iGoldberg et al.l (120021 ) , the 



velocity difference between the two components of this object is expected to be 15-20 km at 
this epoch. Indeed, the separation of the absorption features in our spectrum is 18.5 km s~^ (see 
below), agreeing well with the above estimate. 

Standard data reduction was carried out with the IRAF echelle packag^. Spectra are ex- 
tracted from i ndivid ual sliced images and combined for the slice numbers 2-5 in the definition of 
Tajitsu et al.l (j2012l ) to obtain the best quality spectrum. 



Equivalent widths (W^'s) for isolated absorption lines were measured by fitting Gaussian pro- 
files for each component of the binary. The equivalent widths are used in the following chemical 
abundance analysis by dividing the fraction of the contribution of each component to the continuum 
flux estimated from colors and mass ratio of the two components (§[3]). 



^IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science Foundation. 
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3. Analysis and results 

3.1. Atmospheric Parameters 

In order to determine the atmospheric parameters of the two components of the binary system 
and their contributions to the continuum flux, we adopt the mass ratio of the two components 
{q = m2/mi = 0.89 ifc 0.04), where mi and m2 are the primary and secondary masses, respectively. 



determined by iGoldberg et alj (|2002l ). and calculate the colors (B — V, V — R, and — /) as a 
function of the primary mass using the isochrone ( Demarque et al.ll2004l ) for the age of 12 Gyr. 



The metallicity is assumed to be [Fe/H]= —2.5, which is consistent with the results obtained from 
the following analysis. 

The upper panel of Figure [T] shows the three colors calculated using the isochrone. Solid and 
dotted lines indicate the cases of q = 0.89 and those with 0.04 higher or lower, respectively. Very 
accurate colors of G 166-45 are available, beca use this object is a Landq lt's photometric standard 



star. The colors measured by their recent work (jLandolt &: Uomotd 120071 ) with a reddening correc- 



tion are shown by the horizontal lines in the figure. The reddening correction {E{B — V) = 0.02) 
is estimated from the interstellar Na I D2 line detected in our spectrum (the equivalent width of 
80 mA). 

Comparisons of these colors with the calculations provide constraints on the primary mass of 
this binary system (as well as the secondary mass). The primary mass estimated from the three col- 
ors is O.757±O.OO7M0, in which the error corresponds to the differences of the re sults from the three 



colors . Even if the photometry error (0.015 mag) for G 166-45 estimated by iLandolt &: Uomoto 



(|2007l ) (see the erratum of the paper) and the uncertainty of the mass ratio are included, the error 
size of the primary mass does not change significantly. The other solution that satisfies the colors 
and the mass ratio is the case that the primary is a subgiant (mi ~ 0.82 M0) and the secondary 
is a main-sequence star with similar temperatures (see Fig. [1]). This solution is, however, excluded 
because the luminosity difference of the two components can not be so large, given the strengths 
of the absorption features of the two components. 

The lower panel of Figure [T] shows the effective temperatures of the primary and secondary, 
as a function of the primary mass (assuming the mass ratio), adopted from the isochrone. The 
effective temperature of the primary is well determined to be 6350ibl00 K. The error refiects the 
uncertainty of the mass determination from different color indices and photometric errors, but does 
not include the uncertainty of the temperature scale (depending on model atmospheres). 

By contrast, the uncertainty of the effective temperature of the secondary is relatively large 
(170 K) due to the uncertainty of the mass ratio. We adopt 5830 K for the secondary as the best 
estimate in the following analysis, but also conduct the analyses for Tgfj = 5650 and 6000 K to 
estimate the errors due to the uncertainty of the mass ratio. The surface gravities of the primary 
and the secondary are estimated from the isochrone to be log 5 = 4.4 and 4.6, respectively. 
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The fraction of the contributions of each component to the continuum flux in the i?-band (/: 
/a and /b for the primary and secondary, respectively), in which the Li I resonance line exists, is 
also estimated from the isochrone for the three choices of the secondary's effective temperature. 

The adopted parameter sets (Tcff, logg, and /b) are (5830 K, 4.6, 0.27), (5650 K, 4.6, 0.21) 
and (6000 K, 4.6, 0.34) for the secondary, while Tgfj = 6350 K, log 5 = 4.4, and /a = 1 — /b are 
adopted for the primary (Table 1). 

For comparison purposes, we applied the above procedure to the double-lined sp ectroscopic 
binary CS 22876-032 adop ting the mass ratio of 0.91 (iGonzalez Hernandez et al.l 120081) and colors 



given 



by lNorris et al 



2OOOI ) . The results agree well with those determined by lGonzalez Hernandez et al 
1)20081 ). i.e., 6500 K and 5900 K for the primary and secondary, respectively, with uncertainties of 
about 100 K. 



3.2. Chemical abundances 

Chemical compositions of the two components of this binary system are determined separately 
based on the standard LTE analysis of the equivalent widths obtained by dividing the measured 
equivalent widths from t he spectrum by the /. The ATLAS NEWODF model atmosphere grid 



(jCastelli &: Kuruczll2003l ) assuming a elements excesses is used in the analysis. 



The micro-turbulent velocity (f micro) is determined by demanding no dependence of the abun- 
dance results for individual Fe I lines on their strengths (equivalent widths) . The derived Vmicro is 
1.5 km s~^ for the primary and 1.0 km s~^ for the secondary. 

The results of the abundance analyses are given in Table [T] for eight elements (along with the 
Li abundance determined separately) for the standard case of the stellar parameters. We also give 
the results of Fe abundances for the other parameter sets in the table. The errors given in Table [1] 
include random errors and those due to uncertainties of atmospheric parameters. We estimate 
the random errors in the measurements to be cr/N, where a is the standard deviation of derived 
abundances from individual lines and N is the number of lines used. When only a few lines are 
available, the <j of Fe I is adopted in the estimates. The errors due to the uncertainty of the 
atmospheric parameters {5Te{i= 100 K, 6logg= 0.3, and 5vmicTo= 0.2 km s~^) are also estimated 
and added in quadrature to the random errors. The errors for the secondary could be slightly larger 
if the the uncertainty of effective temperature due to the uncertainty of the mass ratio (~ 170 K) 
is included. 

The Fe abundances derived from individual Fe I lines show no statistically significant depen- 
dence on their lower excitation potentials for the three choices of the stellar parameters. The Fe 
abundance derived from Fe II lines is slightly higher than that from Fe I lines, suggesting the log g 
values adopted in the analysi s are systemati cally too high (by 0.1-0.3 dex), or the Fe I lines suffer 
from non-LTE effects (e.g., Asplundl 2005). We note that the Li line analysis reported in the 
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following subsection is not sensitive to the changes of the adopted logs' and Wmicro- 

The chemical compositions of the two components agree very well with each other. The changes 
of the Fe abundance from Fe I lines by changing the choice of stellar parameter set are only 0.05 dex 
for the primary. The change for the secondary is even smaller, because the effect of the change 
of Tgff on the derived Fe abundance is partially compensated by the change of the contribution of 
the secondary to the continuum flux (adopting a higher Tgg results in a higher Fe abundance in 
general, while a larger /b results in smaller equivalent widths and a lower abundance). 



3.3. Li abundances 

The Li abundances of the two components of G 166-45 are determined by fitting synthetic 
spectra to the observed one. The synthetic spectra are calculated using the same model atmospheres 
used in the above subsection for both components, and are combined taking their /'s into account 
and the velocity shift between them (18.5 km s~^). Comparisons of the synthetic spectra with 
the observed one are depicted in Figure [2 where the doppler correction is made for the primary. 
In the fitting process, the minimum is searched for each component: the wavelength ranges 
6707.6-6708.0 A and 6708.0-6708.4 A are used to determine the Li abundances of the primary 
and secondary, respectively. The 2a error in the fitting is 0.02 and 0.03 dex for the primary and 
secondary, respectively. 

The results for the three choices of mass ratios are given in Table [H Interestingly, the Li 
abundance of the secondary is insensitive to the choice of the mass ratio for the same reason for 
the behavior of the Fe abundance mentioned above. 

The derived Li abundance of the secondary is 0.05-0.15 dex lower than that of the primary. 
In the following discussion on the Li abundance difference between the two components, we take 
the fitting errors (0.02-0.03 dex) and the errors due to the uncertainty of the mass ratio (0.02- 
0.05 dex) into consideration. We note that this does not include the error due to the uncertainty 
of Teff scales, which is as large as 0.1 dex in ^(Li), but is a systematic error affecting the results 
for both components. 



4. Discussion and concluding remarks 



The Li abundance of the primary derived by the above analysis (^(Li)= 2.23ib0.05) well agrees 
with the Spite plateau value, while that of the secondary is slightly lower. Figure [3] (top panel) 
depicts the Li abundances of main-sequence {\.ogg> 4.0: filled symbols) and subgiant {\.ogg< 
4.0: open symbols) stars as a functi on of [Fe/H]. The dat a other than G 166-4 5 are adopted 
from bonzalez Hernandez et al.l (I2OO8I ) for OS 22876-032 and lMelendez et al.l (I2OI0I ) for other stars 
determined by LTE analyses. The possible Li depletion in the main-sequence phase can be inspected 
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from filled symbols shown here. While large scatter of Li abundances is found in [Fe/H]> —2, 
which is due to the lower abundances in stars with lower T^q (see below), there is no distinctive 
scatter in —3 <[Fe/H]< —2 in t his sample. This is not changed even if subgiants (log5f< 4) are 



included. iMelendez et al.l ([20101) suggested the existence of two plateaus with slightly higher and 
lower Li abundances in [Fe/H]> —2.5 and < —2.5, respectively. The Li abundance of G 166-45A 
is apparently in agreement with the lower plateau. This result is, however, not definitive, given a 
possible systematic offset in Li abundance determinations (as large as 0.1 dex) mostly due to Tgff 
scales used. Among the stars in this metallicity range, the Li abundance of G 166-45B is relatively 
low. 

Some scatter of Li abundances is also found in [Fe/H]< —3, though that is smaller than that 
in [Fe/H]> —2. The discrepancy of the two components of CS 22876-032 is consistent with the 
size of the scatter. 

The middle panel of Figure [3] shows the Li abundances as a function of Tefj. A clear depen- 
dence of Li abundances on Tefj is found in Teff< 6000 K. It should be noted that all stars in this 
temperature range, except for CS 22876-032B and two subgiants, are objects with [Fe/H]> —2. 
The difference of Li abundances between objects with Teff~ 6350 K and those with ~5830 K is 
as large as 0.4 dex. The discrepancy of the Li abundances found between the two components of 
G 166-45 is much smaller than this value. This suggests that the depletion of Li in stars with 
Teff~ 5800 K is much smaller at this metallicity ([Fe/H]~ —2.5) than that in [Fe/H]> —2. This 
might be related to the shallower convection zone in sta rs with lower metalli city, however, other 



causes for Li depletion would be required as mentioned by lRichard et al.l (|2005l ). who demonstrated 
that the temperature of the bottom of the convective zone is almost independent of metallicity. It 
should also be noted that subgiants {logg< 4.0) with Tcg< 6000K have Li abundances as high as 
stars with Tcff> 6000 K, hence, the correlation between the Li abundance and T^g discussed above 
is only found in main-sequence stars with high [Fe/H] (> —2). The exception is CS 22876-032B, 
which shows a Li abundance along with the correlation found for stars with [Fe/H]> —2. The rea- 
son for the Li depletion in this object is, however, likely different from the Teff-dependent depletion 
found for [Fe/II]> —2 (see below). 

The bottom panel of Figure [3l sh ows Li abundan c es as a function of stellar mass. Li depletion 



in lower mass stars is discussed by IMelendez et al.l (|2010l ). Excluding the stars with [Fe/H]> 



-2 (indicated by triangles), which show a tight correlation between Li ab undances and T p ff, the 



dependence of the Li abundance on stellar mass is not very clear. Indeed, IMelendez et al.l (|2010l ) 
reported that the slopes of Li abundances as a function of stellar masses is significant at only 
la level in [Fe/H]< —2.5, which increases to the 3a level by including the binary CS 22876-032. 
However, if stars with lowest metallicity ([Fe/II]< —3) are selected, a correlation seems to appear 
again. In order to demonstrate this, stars in this metallicity range are indicated by over-plotting 
open circles (except for CS 22876-032) in the bottom panel of the figure. 



On the other hand, no clear correlation is found between the Li abundances and stellar masses 
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for objects with —3 <[Fe/H]< —2 (squares without open circles). Our study of the Li abundance 
of G 166-45 B extends this trend to a mass below 0.7 Mq. 

The dependence of Li abundances on stellar mass and metallicity is nonlinear. In [Fe/H]> —2, 

Li abundances arc lower in main-sequence stars with lower T^s (lower mass) in T^g^ 6000 K 
(M < O.SMq). a similar dependence is also found in [Fe/H]< —3, possibly in stars with even 
higher T^s (< 6300 K, which corresponds to M < O.SMq). On the other hand, such a dependence 
is not clearly seen in —3 <[Fe/H]< —2. The double-lined spectroscopic binaries CS 22876-032 and 
G 166-45 most clearly represent these trends in [Fe/H]< —3 and —3 <[Fe/H]< —2, respectively. 
Since the initial abundances of the two components of a binary system can be assumed to be 
the same, the discrepancy of the Li abundances found between the two components indicates a 
depletion of Li in the secondary. This suggests that the scatter of Li found in the lowest metallicity 
range is due to depletion of Li in some object, which is not clearly seen in —3 <[Fe/H]< —2. 

Facilities: Subaru(HDS). 

We are grateful to M. Parthasarathy for his useful comments for the manuscript. W.A. was 
supported by the Grants-in-Aid for Science Research of JSPS (20244035). 
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Fig. 1. — Upper panel: Colors iy — I,B — V,V — R) calculated from the isochrones assuming mass 
ratios of g = 0.89 ± 0.04. The reddening corrected colors of G 166-45 are shown by horizontal 
lines. The vertical solid and dashed lines indicate the adopted primary mass and its possible range 
estimated from the three colors. Lower panel: effective temperatures of primary (thick solid line) 
and secondary (thin solid line and dotted lines) as a function of assumed primary mass adopting 
the above mass ratios. 
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Table 1. STELLAR PARAMETERS AND CHEMICAL ABUNDANCES 
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Ba II 


-0.47 


-0.20 


1 


0.14 










Li I 


2.23 




1 


0.02^ 


2.11 




1 


0.03^^ 


q = 0.93 


reff=6350 


/a = 0.66 








=6000 K, /b = 0.34 






Fe I 


5.11 


-2.39 


35 


0.10 


5.09 


-2.42 


39 


0.11 


Fe II 


5.13 


-2.37 


3 


0.14 


5.19 


-2.31 


2 


0.15 


Li I 


2.27 




1 


0.02^ 


2.12 




1 


0.03^ 


q = 0.85 


Tefr=6350 


/a = 0.79 






TeS= 


=5670 K, /b = 0.21 






Fc I 


5.00 


-2.50 


35 


0.10 


5.14 


-2.36 


39 


0.11 


Fe II 


5.04 


-2.46 


3 


0.14 


5.43 


-2.07 


2 


0.15 


Li I 


2.18 




1 


0.02^ 


2.13 




1 


0.03^ 



^The errors of Li abundances represent random errors in the fitting of synthetic 
spectra. 



